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Vector theories with spontaneous Lorentz violation, known as bumblebee models, are examined in 
flat spacetime using a Hamiltonian constraint analysis. In some of these models, Nambu-Goldstone 
modes appear with properties similar to photons in electromagnetism. However, depending on 
the form of the theory, additional modes and constraints can appear that have no counterparts in 
electromagnetism. An examination of these constraints and additional degrees of freedom, including 
their nonlinear effects, is made for a variety of models with different kinetic and potential terms, and 
the results are compared with electromagnetism. The Hamiltonian constraint analysis also permits 
an investigation of the stability of these models. For certain bumblebee theories with a timelike 
vector, suitable restrictions of the initial-value solutions are identified that yield ghost-free models 
with a positive Hamiltonian. In each case, the restricted phase space is found to match that of 
electromagnetism in a nonlinear gauge. 



I. INTRODUCTION 

Investigations of quantum-gravity theories have un- 
covered a variety of possible mechanisms that can lead 
to Lorentz violation. Of these, the idea that Lorentz 
symmetry might be spontaneously broken [l| is one of 
the more elegant. Spontaneous Lorentz violation occurs 
when a vector or tensor field acquires a nonzero vacuum 
expectation value. The presence of these background val- 
ues provides signatures of Lorentz violation that can be 
probed experimentally. The theoretical framework for 
their investigation is given by the Standard-Model Exten- 
sion (SME) Experimental searches for low-energy 
signals of Lorentz violation have opened up a promising 
avenue of research in investigations of quantum-gravity 
phenomenology 0, [|. 

Theories with spontaneous Lorentz violation can also 
exhibit a variety of physical effects due to the appear- 
ance of both Nambu-Goldstone (NG) and massive Higgs 
modes 0, 0, [| • In the context of a gravitational theory, 
these effects include modifications of gravitational prop- 
agation, as well as altered forms of the static Newtonian 
potential, both of which may be of interest in theoreti- 
cal investigations of dark energy and dark matter. Many 
investigations to date have concentrated on the case of 
a vector field acquiring a nonzero vacuum value. These 
theories, called bumblebee models [jj, [9|, tlQj] , are the sim- 
plest examples of field theories with spontaneous Lorentz 
breaking. Bumblebee models can be defined with differ- 
ent forms of the potential and kinetic terms for the vector 
field, and with different couplings to matter and gravity 
[H E3, 0, El OS El EE m- The Y can be considered 
as well in different spacetime geometries, including Rie- 
mann, Riemann-Cartan, or Minkowski spacetimes. 

Much of the interest in bumblebee models stems from 
the fact that they are theories without local ?7(1) gauge 
symmetry, but which nonetheless allow for the propaga- 
tion of massless vector modes. Indeed, one idea is that 
bumblebee models, with appropriate kinetic and poten- 
tial terms, might provide alternative descriptions of pho- 



tons besides that given by local U(l) gauge theory. In 
this scenario, massless photon modes arise as NG modes 
when Lorentz violation is spontaneously broken. How- 
ever, in addition to lacking local U(l) gauge invariance, 
bumblebee models differ from electromagnetism (in flat 
or curved spacetime) in a number of other ways. For ex- 
ample, the kinetic terms need not have a Maxwell form. 
Instead, a generalized form as considered, for example, in 
vector-tensor theories of gravity can be used, though typ- 
ically this may involve the introduction of ghost modes 
into the theory. Further differences arise due to the pres- 
ence of a potential term V in the Lagrangian density for 
bumblebee models. It is this term that induces sponta- 
neous Lorentz breaking. It can take a variety of forms, 
which may involve additional excitations due to the pres- 
ence of massive modes or Lagrange-multiplier fields that 
have no counterparts in electromagnetism. 



The goal of this paper is to investigate further the ex- 
tent to which bumblebee models can be considered as 
equivalent to electromagnetism or as containing electro- 
magnetism as a subset theory. This question is examined 
here in flat spacetime. While gravitational effects are 
a feature of primary interest in bumblebee models, any 
equivalence or match to electrodynamics would presum- 
ably hold as well in an appropriate flat-spacetime limit. 
In a Minkowski spacetime, the main differences between 
bumblebee models and electromagnetism are due to the 
nature of the constraints imposed on the field variables 
and in the number of physical degrees of freedom per- 
mitted by the theory. To investigate these quantities, a 
Hamiltonian constraint analysis [19l . I20L HH, Hll is used. 
This approach is particularly well suited for identifying 
the physical degrees of freedom in a theory with con- 
straints. It can be carried out exactly with all nonlinear 
terms included. It also permits examination of the ques- 
tion of whether the Hamiltonian is bounded from below 
over the constrained phase space. 
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II. BUMBLEBEE MODELS AND 
ELECTROMAGNETISM 

Bumblebee models are field theories with spontaneous 
Lorentz violation in which a vector field acquires a 
nonzero vacuum value. For the case of a bumblebee 
field Bp coupled to gravity and matter, with general- 
ized quadratic kinetic terms involving up to second-order 
derivatives in B^, and with an Einstein-Hilbert term for 
the pure-gravity sector, the Lagrangian density is given 
as 

lOTTLr 

+ ^t 3 D^D„B» -V(B^Tb 2 )+£M- (1) 

In this expression, b 2 > is a constant, and in Riemann 
spacetime B^ u = d^B v ~d v B^. The quantities <7i, 02, ri, 
T2, and T3 are fixed constants that determine the form of 
the kinetic terms for the bumblebee field. The term Cm 
represents possible interaction terms with matter fields 
or external currents. The potential V(B fl B^ 1 ^fb 2 ) has a 
minimum with respect to its argument or is constrained 
to zero when 

B^ T b 2 = 0. (2) 

This condition is satisfied when the vector field has a 
nonzero vacuum value 

B^ = (BJ = b„ (3) 

with b^ = ±b 2 . It is this vacuum value that sponta- 
neously breaks Lorentz invariance. 

There are many forms that can be considered for the 
potential V(B fl B ,J ' =p b 2 ). These include functionals in- 
volving Lagrange- multiplier fields, as well as both poly- 
nomial and nonpolynomial functionals in {B^B^ =p b 2 ) 
[H> [Hi- I n this work, three limiting-case examples are 
considered. They represent the dominant leading-order 
terms that would arise in an expansion of a general scalar 
potential V, comprised of vector fields B^ 7 which are 
not simply mass terms. They include examples that 
are widely used in the literature. The first introduces 
a Lagrange-multiplier field A and has a linear form, 

V = \(B^ T b 2 ). (4) 

which leads to the constraint ^ appearing as an equa- 
tion of motion. The second is a smooth quadratic poten- 
tial 

V = \n{B^Tb 2 )\ (5) 

where n is a constant. The third again involves a 
Lagrange-multiplier field A, but has a quadratic form, 

V = ±\(B^ T b 2 ) 2 . (6) 



With this form, the Lagrange multiplier field A decouples 
from the equations of motion for B^. 

The model given in ([I]) involving a vacuum- valued vec- 
tor has a number of features considered previously in the 
literature. For example, with the potential V and the 
cosmological constant A excluded, the resulting model 
has the form of a vector-tensor theory of gravity consid- 
ered by Will and Nordvedt [HHl]. Models with poten- 
tials (@| and ([5]) inducing spontaneous symmetry break- 
ing were investigated by Kostelecky and Samuel (KS) 
while the potential ([6]) was recently examined in Q • The 
special cases with a nonzero potential V, n = 1, and 
f 1 = 02 = T2 = T3 = are the original KS bumblebee 
models [l[ . Models with a linear Lagrange-multiplier po- 
tential (|4]), <7i = (T2 = 0, but arbitrary coefficients n, T2, 
and T3 are special cases (with a fourth-order term in B^ 
omitted) of the models described in Ref. 12]. 

Since bumblebee models spontaneously break Lorentz 
and diffeomorphism symmetry, it is expected that mass- 
less Nambu-Goldstone (NG) and massive Higgs modes 
should appear in these theories. The fate of these modes 
was recently investigated in @] • The example of a KS 
bumblebee was considered in detail. It was found that 
for all three potentials (J3J, ©, and ([6]), massless NG 
modes can propagate and behave essentially as photons. 
However, in addition, it was found that massive modes 
can appear that act as additional sources of energy and 
charge density. In a linearized and static limit of the KS 
bumblebee, it was shown that both the Newtonian and 
Coulomb potentials for a point particle are altered by the 
presence of a massive mode. Nonetheless, with suitable 
choices of initial values, which limit the phase space of the 
theory, solutions equivalent to those in Einstein-Maxwell 
theory can be obtained for the KS bumblebee models. 

Bumblebee models with other (non-Maxwell) values of 
the coefficients t±, T2, and T3 are expected to contain 
massless NG modes as well. However, in this case, since 
the kinetic terms arc different, a match with electrody- 
namics is not expected. The non-Maxwell kinetic terms 
alter the constraint structure of the theory significantly, 
and a different number of physical degrees of freedom can 
emerge. 

To compare the constraint structures of different types 
of bumblebee models with each other and with electro- 
dynamics, the flat-spacetime limit of (fT|) is considered. 
The Lagragian density in this case reduces to 

C = —nB^B^ + ^d„B v a^B v 

+ \r 3 d^d u B" - V(B^ ± b 2 ) - B^. (7) 

For simplicity, interactions consisting of couplings with 
an externally prescribed current J M are assumed, and 
a Minkowski metric rj^ in Cartesian coordinates with 
signature (+, — , — , — ) is used. 

Following a Lagrangian approach, second-order differ- 
ential equations of motion for B^ are obtained. They 
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(n + r 3 ) [ps„ - d^d v B v ] 

-{T 2 + T 3 )nB, l -2V'B f ,- J M = 0. (8) 

Here, V denotes variation of the potential V(X) with 
respect to its argument X. Since the NG modes stay in 
the minimum of the potential, a nonzero value of V indi- 
cates the presence of a massive-mode excitation. Taking 
the divergence of these equations gives 



[(r 2 + t 3 )UB„. + 2V'B I1 + J„] = 0. 



(9) 



Clearly, as expected, with V = V = 0, t\ = 1, and the 
remaining coefficients set to zero, the equations of motion 
reduce to those of electrodynamics, and reduces to the 
statement of current conservation. However, if a nonzero 
potential with V ^ 0, or if arbitrary values of n, r%, T3 
are allowed, then a modified set of equations holds. 

In flat spacetime, the KS bumblebee has a nonzero 
potential V and coefficients T\ = 1, and r 2 = T3 = 0. Its 
equations of motion evidently have a close resemblance 
to those of electrodynamics. The main difference is that 
the KS bumblebee field itself acts nonlinearly as a source 
of current. Equation © shows that the matter current 
J p combines with the term 2V'B /1 to form a conserved 
current. 

Interestingly, if the matter current J M is set to zero, 
and a linear Lagrange-multiplier potential @ is used, 
the KS model in flat spacetime reduces to a theory con- 
sidered by Dirac long before the notion of spontaneous 
symmetry breaking had been introduced [25]. Dirac in- 
vestigated a vector theory with a nonlinear constraint 
identical to ([2]) with the idea of finding an alternative ex- 
planation of electric charge. In his model, gauge invari- 
ance is destroyed, and conserved charge currents appear 
only as a result of the nonlinear term involving V' for the 
Lagrange-multiplier potential. Dirac did not, however, 
propose a theory of Lorentz violation. A vacuum value 
bft was never introduced, and with J M = no Lorentz- 
violating interactions with matter enter in the theory. 

The idea that the photon could emerge as NG modes in 
a theory with spontaneous Lorentz violation came more 
than ten years after the work of Dirac. First, Bjorken 
proposed a model in which collective excitations of a 
fermion field could lead to composite photons emerging as 
NG modes [26] . The observable behavior of the photon in 
this original model was claimed to be equivalent to elec- 
trodynamics. Subsequently, Nambu recognized that the 
constraint ||2J| imposed on a vector field could also lead 
to the appearance of NG modes that behave like photons 
[27l ] . He introduced a vector model that did not involve a 
symmetry-breaking potential V. Instead, the constraint 
@ was imposed as a nonlinear U(l) gauge- fixing condi- 
tion directly at the level of the Lagrangian. The resulting 
gauge-fixed theory thus contained only three indepen- 
dent vector-field components in the Lagrangian. Nambu 
demonstrated that his model was equivalent to electro- 
magnetism and stated that the vacuum vector can be 
allowed to vanish to restore full Lorentz invariance. 



In contrast to these early models, the KS bumblebee 
was proposed as a theory with physical Lorentz viola- 
tion. Even if the NG modes are interpreted as photons 
in the KS model, and no massive modes are present, in- 
teractions between the vacuum vector & M and the matter 
current provide clear observable signals of physical 
Lorentz violation. However, the presence of a potential 
V also allows additional degrees of freedom to enter in 
the KS model. If arbitrary values of the coefficients n, 
T2, and T3 are permitted as well, the resulting theory can 
differ substantially from electromagnetism. 

Since many of these models contain unphysical modes, 
either as auxiliary or Lagrange-multiplier fields, con- 
straint equations are expected to hold. It is the nature of 
these constraints that determines ultimately how many 
physical degrees of freedom occur in a given model. With 
Dirac's Hamiltonian constraint analysis, a direct proce- 
dure exists for determining the constraint structure and 
the number of physical degrees of freedom in these mod- 
els. 



III. HAMILTONIAN CONSTRAINT ANALYSIS 

Given a Lagrangian density C describing a vector field 
B^, the canonical Hamiltonian density is H — H^doB^ — 
C, where the canonical momenta are defined as 



SC 



5{d B, 



(10) 



If additional fields, e.g., Lagrange multipliers A, are con- 
tained in the theory, additional canonical momenta for 
these quantities are defined as well, e.g., n( A ). (Note: 
here A is not a spacetime index) . In the Hamiltonian ap- 
proach, time derivatives of a quantity / are computed by 
taking the Poisson bracket with the Hamiltonian H, 



(11) 



The second term is needed with quantities that have ex- 
plicit time dependence, e.g., an external current J M . 

In Dirac's constraint analysis, primary and secondary 
constraints are determined, and these are identified as ei- 
ther first-class or second-class. In the phase space away 
from the constraint surface, the canonical Hamiltonian is 
ambiguous up to additional multiples of the constraints. 
An extended Hamiltonian is formed that includes multi- 
ples of the constraints with coefficients that can be de- 
termined, or in the case of first-class constraints, remain 
arbitrary. It is the extended Hamiltonian that is then 
used in (jlip to determine the equations of motion for the 
fields and conjugate momenta. 

A system of constraints is said to be regular if the Ja- 
cobian matrix formed from variations of the constraints 
with respect to the set of field variables and conjugate 
momenta has maximal rank. If it does not, the system 
is said to be irregular, and some of the constraints are 
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typically redundant. Dirac argued that theories with pri- 
mary first-class constraints have arbitrary or unphysical 
degrees of freedom, such as gauge degrees of freedom. 
These types of constraints therefore allow removal of two 
field or momentum components. Dirac conjectured that 
this is true as well for secondary first-class constraints. 
Based on this, a counting argument can be made. It 
states that in a theory with n field and n conjugate- 
momentum components, if there are ri\ first-class con- 
straints and n 2 second-class constraints, the number of 
physical independent degrees of freedom is n — m — n 2 /2. 
(Note: it can be shown that n 2 is even). This counting 
argument based on Dirac's conjecture holds up well for 
theories with regular systems of constraints. However, 
counterexamples are known for irregular systems [22] . 

Once the unphysical modes have been eliminated, by 
applying the constraints and/or imposing gauge condi- 
tions, the evolution of a physical system is determined 
by the equations of motion for the physical fields and 
momenta, subject to initial conditions for these quanti- 
ties. Any bumblebee theory that has additional degrees 
of freedom in comparison to electrodynamics must there- 
fore specify additional initial values. The subsequent evo- 
lution of the extra degrees of freedom typically leads to 
effects that do not occur in electrodynamics. However, 
in some cases, equivalence with electrodynamics can hold 
in a subspace of the phase space of the modified theory. 
For this to occur, initial values must exist that confine 
the evolution of the theory to a region of phase space that 
matches electrodynamics in a particular choice of gauge. 

In general, the stability of a theory, e.g., whether the 
Hamiltonian is positive, depends on the initial values and 
allowed evolution of the physical degrees of freedom. As 
discussed in the subsequent sections, most bumblebee 
models contain regions of phase space that do not have a 
positive definite Hamiltonian, though in some cases, re- 
stricted subspaces can be found that do maintain TL > 0. 
In a quantum theory, instability in any region of the clas- 
sical phase space might be expected to destabilize the 
full theory. However, bumblebee models, with gravity 
included, are intended as effective theories presumably 
emerging at or below the Planck scale from a more fun- 
damental (and unknown) quantum theory of gravity. In 
this context, quantum-gravity effects might impose ad- 
ditional constraints leading to stability. However, in the 
absence of a fundamental theory, the question of the ulti- 
mate stability of bumblebee models cannot be addressed. 
For this reason, in the subsequent sections, only the be- 
havior of bumblebee models in classical phase space is 
considered. 

The following sections apply Dirac's constraint analy- 
sis to a number of different bumblebee models, including 
the KS bumblebee as well as more general cases with ar- 
bitrary values of the coefficients t\, t 2 , T3. Since much of 
the literature has focused on the case of a timelike vec- 
tor B^, this restriction is assumed throughout this work 
as well. With this assumption, there always exists an 
observer frame in which rotational invariance is main- 



tained and only Lorentz boosts are spontaneously bro- 
ken. For each type of model to be considered, all three 
of the potentials in ©, and ^ are considered. For 
comparison (and use as benchmarks), electromagnetism 
and the theory of Nambu are considered as well. In each 
case, the explicit form of the Lagrangian is obtained from 
(|7|) by inserting appropriate values for V, t\, t 2 , and T3, 
and the conjugate momenta and Hamiltonian are then 
computed. For example, electrodynamics is obtained by 
setting V = 0, t\ = 1, and t 2 = T3 = 0. Conventional no- 
tation sets = A^ and _B M „ = . The Hamiltonian is 
given in terms of the four fields and their conjugate 
momenta n M . The Lagrangian in Nambu's model also 
starts with these same values (allowing U(l) invariance). 
However, in this case, one component of A^ is eliminated 
in terms of the remaining three, using the nonlinear con- 
dition in ([2]). For the case of a timelike vector, the sub- 
stitution A = (b 2 + A 2 ) 1 / 2 is made directly in C. The re- 
sulting Hamiltonian in Nambu's model therefore depends 
only on three fields Aj and three conjugate momenta 
IP. In contrast, bumblebee models are defined with a 
nonzero potential V and have Hamiltonians that depend 
on all four fields B^ and their corresponding conjugate 
momenta n M . Examples with a Lagrange- multiplier po- 
tential involve a fifth field A and its conjugate momentum 
H( x \ However, in examples with a smooth quadratic po- 
tential, there is no Lagrange multiplier, and the relevant 
fields and momenta are B^ and n M . 



A. Electromagnetism 

The conjugate momenta in electrodynamics are 

= OoAj - djA , n° = 0. (12) 

The latter constitutes a primary constraint, 4>\ — n° w 0. 
It leads to a secondary constraint, <f> 2 = cjj-EP — J° ~ 0, 
which is Gauss' law, since IP can be identified as the elec- 
tric field components and J° is the charge density. In 
these expressions and below, Dirac's weak equality sym- 
bol is used to denote equality on the submanifold de- 
fined by the constraints [22|. Both of the constraints <fii 
and <p2 are first-class, indicating that there are gauge or 
unphysical degrees of freedom. Following Dirac's count- 
ing argument, there should be n—ri\— n 2 /2 = 4—2—0 = 2 
independent physical degrees of freedom. These are the 
two massless transverse photon modes. 

The canonical Hamiltonian in electrodynamics is 

TL = \{W) 2 + WdjAo + ±(F jk ) 2 + A^. (13) 

In the presence of a static charge distribution, with 
J p = (p, J) = (p(x),0), no work is done by the exter- 
nal current, and the Hamiltonian is positive definite. To 
observe this, integrate by parts and use the constraint <p 2 
(Gauss' law) to show that TL = ±(n 1 ) 2 + \{F jk ) 2 > 0. 



5 



The equations of motion for the fields and momenta 
11^ obtained from the extended Hamiltonian contain ar- 
bitrary functions due to the existence of the first-class 
constraints. These can be eliminated by imposing gauge- 
fixing conditions. The evolution of the physical degrees 
of freedom, subject to a given set of initial values, is then 
determined for all time. 



B. Nambu's Model 

The starting point for Nambu's model [13] is the con- 
ventional Maxwell Lagrangian with U(l) gauge invari- 
ance and a conserved current J M . For the case of a time- 
like vector Ap, the condition Aq = (b 2 + A 2 ) 1 / 2 is sub- 
stituted directly into the Lagrangian as a gauge-fixing 
condition. The result is 



C = ^(doA,) 2 + l(d 3 (b 2 + A 2 )^ 2 ) 2 

md,A k )(d k A,) - (d 3 (b 2 

-{b 2 + Al) 1 ' 2 J° -AiJi. 



Alf' 2 



(14) 



Nambu claimed that this theory is equivalent to electro- 
magnetism in a nonlinear gauge. He argued that a U(l) 
gauge transformation exists that transforms an electro- 
magnetic field in a standard gauge into the field A^ obey- 
ing the nonlinear gauge condition A^A^ = b 2 . 
The Hamiltonian in Nambu's model is 



H 



\{W) 2 + \{F ]k ) 2 + Wd, J {b 2 + Alf' 2 
+ (b 2 + A 2 ) 1 / 2 J° 



It depends on three field components Aj 
jugate momenta IF' = OqAj — dj(b 2 



(15) 

and their con- 
Al) 1 / 2 . In this the- 
ory, there are no constraints, and therefore application of 
Dirac's counting argument says that there arc three phys- 
ical degrees of freedom, which is one more than in elec- 
tromagnetism. An extra degree of freedom arises because 
gauge fixing at the level of the Lagrangian causes Gauss' 
law, djlP —J = 0, to disappear as a constraint equation. 
A smilar disappearance of Gauss' law is known to occur 
in electrodynamics in temporal gauge (with Aq = sub- 
stituted in the Lagrangian) 28]. Indeed, the linearized 
limit of Nambu's model with a timelike vector field is 
electrodynamics in temporal gauge. 

Observe that with with J = and using integration 
by parts, the Hamiltonian can be rewritten as 



n = f + -(F jk ) 2 - (d 3 W - J°)(b 2 + A 2 ). (16) 



In the absence of a constraint enforcing Gauss' law, Tt 
need not be positive definite. For example, if the ex- 
tra degree of freedom in Aj causes large deviations from 
Gauss' law, which are not forbidden by any constraint, 
then negative values of Ti. can occur. 

However, equivalence between Nambu's model and 
electrodynamics can be established by restricting the 



phase space in Nambu's theory. To see that this follows, 
consider the equations of motion in Nambu's model, 



A j =rt+d j (b 2 + A 2 )V 2 , 



IF = d k d k A } - d 3 d k A k - <9,n 



(17) 



A' 



(b 2 + A 2 y/ 2 

P. (18) 



(b 2 + A 2 y/ 2 

Taking the spatial divergence of (fl~8|) and using current 
conservation yields the nonlinear relation 



d (djU j -J°) = -dj 



{din 1 - j°) 



(b 2 + A 2 y/ 2 



(19) 



This equation shows that if Gauss' law, {djW — J°) = 0, 
holds at t = 0, then d (djlp - J°) = as well at t = 0. 
Together these conditions and Eq. (fl9|) are sufficient to 
show that Gauss' law then holds for all time. From this it 
follows that H is positive over the restricted phase space, 
which matches that of electrodynamics in a nonlinear 
gauge. Thus, by restricting the phase space to solutions 
with initial values obeying Gauss' law, the equivalence of 
Nambu's model with electromagnetism is restored. 



C. KS Bumblebee Model 

KS bumblebee models 1] in flat spacetime have a 
Maxwell kinetic term and a nonzero potential V . The 
choice of a Maxwell form for the kinetic term is made to 
prevent propagation of the longitudinal mode of as a 
ghost mode. The KS Lagrangian is obtained from (0 by 
setting t\ = 1 and t-i = T3 = 0. The constraint structures 
for models with each of the three potentials (|4]) - ([6]) are 
considered. For definiteness, the case of a timelike vector 
is assumed. 



1. Linear Lagrange-Multiplier Potential 

With a linear Lagrange- multiplier potential ((4]) , an ad- 
ditional field component A is introduced in addition to 
the four fields Bq and Bj. The conjugate momenta are 

n° = n( A > - 0, w = d Bi - d.Bo, (20) 

and the canonical Hamiltonian is 
H = i(ff) 2 + Wd t B Q + ±(d^) 2 - ^(djB^Bj) 

+\{B 2 -B 2 -b 2 ) + B^. (21) 
Four constraints are identified as 

0! = n° (22) 

02 = n( A ) (23) 

03 = d i n i - 2XB - J° (24) 

(25) 



(B 2 -B 2 - b 2 ) 
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The constraints 0i and 02 are primary, while 03 and 04 
are secondary. All four are second-class. 

Applying Dirac's algorithm to determine the number 
of independent degrees of freedom gives n — n\ — 77,2/2 = 
5 — 2 — 2/2 = 3. Hence, there is an extra degree of freedom 
in the KS bumblebee model in comparison to electrody- 
namics. It arises due the presence of the extra field A and 
the changes in the types of constraints. Unlike electro- 
magnetism, there are no first-class constraints in the KS 
bumblebee, which reflects the lack of gauge invariance. 
The constraint 03 gives a modified form of Gauss' law in 
which the combination 2XBq acts as a source of charge 
density. Since V = A in this example, any excitation 
of the field A is away from the potential minimum and 
therefore acts effectively as a massive Higgs mode Q- 
In curved spacetime, such a mode can modify both the 
gravitational and electromagnetic potentials of a point 
particle. However, here, in flat spacetime, the presence 
of A leads only to modifications of the Coulomb potential. 

The Hamiltonian with J = reduces, after using 03, 
04, and integration by parts, to 



l -{B jk ) 2 -2\B 2 . (26) 



The full phase space of the theory on the constraint sur- 
face includes regions in which TL is negative due to the 
presence of the additional degree of freedom. For ex- 
ample, consider the case with J° — and initial val- 
ues [11] Bj = d <p{x) and IP = -d 3 {b 2 + (<9 fc 0) 2 ) 1/2 
at t = 0, where 4>(x) is an arbitrary time-independent 
scalar. These give B jk = and B = (b 2 + (<9,0) 2 ) 1/2 at 
t = 0. Inserting these initial values in (f2"6"]l reduces the 
Hamiltonian to TL — — MlP) 2 at t = 0. The correspond- 
ing initial value for A is 

A = -\(b 2 + (9,0) 2 )- 1 / 2 [V 2 (6 2 + {d^) 2 ) 1 ' 2 ] . (27) 

Evidently, the Hamiltonian in the classical KS bumble- 
bee model can be negative when nonzero values of A are 
allowed. 

However, if initial values are chosen that restrict the 
phase space to values with A = 0, the resulting solutions 
for the vector field and conjugate momentum are equiv- 
alent to those in electromagnetism in a nonlinear gauge. 
Examination of the equation of motion for A, 

^^(A^-^-^-A^^+^o), 

(28) 

reveals that if the current J M is conserved, and A = 
at time zero, then A will remain zero for all time. The 
Hamiltonian in this case is positive. The equations of 
motion for Bj and W are 

Bj = W + djB , (29) 
fP = d k d k Bj - djd k B k - jl+2\Bj (30) 

With A — 0, these combine to give the usual Maxwell 
equations describing massless transverse photons. The 



third component in Bj is an auxiliary field that is con- 
strained by the usual form of Gauss' law when A = 0. 
Note, however, that even with the phase space restricted 
to regions with A = 0, the matter sector of the theory 
will exhibit signatures of the spontaneous Lorentz viola- 
tion through the interaction of the vacuum value 6 M with 
the matter current J M . 

It is clear from these results, that conservation of the 
matter current J M is necessary for the stability of the 
KS bumblebee model. Note, however, that the theories 
lack local U(l) gauge invariance and that the current 
conservation could arise simply from matter couplings 
that are invariant under a global U(l) symmetry. As a 
result, photons in the KS bumblebee model appearing as 
NG modes are due to spontaneous Lorentz breaking, not 
local U(l) gauge invariance. For further discussion of the 
bumblebee currents, including in the presence of gravity, 
see Ref. Q • In that work, there is also further discussion 
of the fact that the Lagrange-multiplier field can act as a 
source of charge density in the KS bumblebee model and 
that there can exist solutions (with nonzero values of A) 
in which the field lines converge or become singular, even 
in the absence of matter charge. This behavior has been 
referred to in the literature as the formation of caustics in 
the KS model. However, as described in Q, it is simply a 
natural consequence of the fact that the bumblebee fields 
themselves act as sources of current. Moreover, with the 
phase space restricted to regions with A = 0, the only 
singularities appearing for the case of a timelike vector 
B^ are those due to the presence of matter charge as in 
ordinary electrodynamics with a 1/r potential. 



2. Quadratic Smooth Potential 

A similar analysis can be performed for a KS bumble- 
bee with the smooth quadratic potential defined in (|5|). 
The parameter k appearing in V is a constant. There- 
fore, in this case, there are four fields Bo, Bj, and their 



four conjugate momenta, 



n° = 0, = doBj - djB . 



There are two constraints, 

01 = H° 

02 = djW - 2kB (B { 



B 2 - b 2 ) - J°, 



(31) 



(32) 
(33) 



where 0i is primary, 02 is secondary, and both are 
second-class. Dirac's counting argument says there are 
n — n\ — 77.2/2 = 4 — — 2/2 = 3 independent degrees 
of freedom, which again is one more than in electromag- 
netism. 

The condition §2$ does not occur as a constraint in this 
case. Instead, an extra degree of freedom appears as a 
massive Higgs excitation V' = 2kB (Bq - B 2 -b 2 ) 
away from the potential minimum. The constraint 02 
yields a modified version of Gauss' law, showing that the 
massive mode acts as a source of charge density. 
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The stability of the Hamiltonian with J = can be 
examined. Using the constraints and integration by parts 
gives 

K = ^(ff') 2 +i(5 ife ) 2 -i K (3B 2 +S?+6 2 )(5 2 -B^-6 2 ), 

(34) 

which evidently is not positive over the full phase space. 
If a nonzero massive mode proportional to (Bq — B| — 6 2 ) 
is present, negative values of Ji can occur. 

However, equivalence to electrodynamics does hold in 
a restricted region of phase space. To verify this, consider 
the equations of motion, 

2kB = (3B 2 - 5? - & 2 )" 1 [4 K B B k (Tl k +d k B ) 

+2nd k [B k (B 2 - Bf -b 2 )}+ d^] , (35) 
Bj = EP + 8jB , (36) 
IT = d 3 W - J° - 2kB q {BI -B] -b 2 ), (37) 
I I = d k d k B, - iUh lh, 

+2nB 3 (B 2 -B 2 k -b 2 )-J j . (38) 

Combining these gives 

nd (B 2 - B) - b 2 ) = (3B 2 - B 2 - b 2 )- 1 

x [2 K B Q d k [B k (B 2 - Bf - b 2 )] + B d^ 
-2k{B 2 -B 2 - 6 2 )B i (n i + diB )] . (39) 

This equation reveals that if the current J^ 1 is conserved 
and (B 2 - B 2 - b 2 ) = at t = 0, then (B% - B 2 - b 2 ) = 
for all time. Therefore, with these conditions imposed, 
the massive mode never appears, the Hamiltonian is pos- 
itive, and the phase space is restricted to solutions in 
electromagnctism in the nonlinear gauge @. 

In theories with a nonzero massive mode, the size of the 
mass scale nb 2 becomes relevant. For very large values, 
perturbative excitations that go up the potential min- 
imum would be expected to be suppressed. Since the 
mass scale associated with spontaneous Lorentz violation 
is presumably the Planck scale, its appearance necessar- 
ily brings gravity into the discussion. It is at the Planck 
scale where quantum-gravity effects might impose addi- 
tional constraints that could maintain the overall stabil- 
ity of the theory. At sub-Planck energies, massive-mode 
excitations have been shown to exert effects on classical 
gravity. For example, as shown in Ref. [7|, the gravita- 
tional potential of a point particle is modified. However, 
in the limit where the mass of the massive mode becomes 
exceptionally large, it was found for the case of the KS 
bumblebee model that both the usual Newtonian and 
Coulomb potentials are recovered. 

3. Quadratic Lagrange- Multiplier Potential 

The KS bumblebee model with a quadratic Lagrange- 
multiplier potential ^ involves five fields A and B^. In 



a Lagrangian approach, the constraint ([2]) follows from 
the equation of motion for A. The on-shell equations of 
motion for B^ are the same as in electromagnctism. In 
this case, the field A decouples and does not act as a 
source of charge density. On shell, the potential obeys 
V' = 0, current conservation 9 M J M = holds, and there 
is no massive mode. This model provides an example of a 
theory with physical Lorentz violation due to the matter 
couplings with J M . Nonetheless, in the electromagnetic 
sector, the theory is equivalent to electromagnctism in 
the nonlinear gauge 

However, the Hamiltonian formulation of this model 
involves an irregular system of constraints |22j. Thus, 
depending on how the constraints are handled, Dirac's 
counting algorithm might not apply and equivalence with 
the Lagrangian approach may not hold. The conjugate 



momenta are 

n° = 0, (40) 

W = d Bj - djB , (41) 

n (A) = 0. (42) 

From these, four constraints can be identified, 

fa = n°, (43) 

4>2 = rt x \ (44) 

03 = djW - 2XB (B 2 - B) - b 2 ) - J°, (45) 

fa = -\(Bl-B 2 -b 2 f. (46) 



With 04 w 0, the constraint surface is limited to fields 
obeying (B 2 — B 2 — b 2 ) = 0, and 03 reduces to Gauss' 
law. In this case, 02 and 04 can be identified as first- 
class, while 0i and 3 are second-class. Dirac's counting 
algorithm then states that there are n — n\ — 112/ 2 = 
5 — 2 — 2/2 = 2 independent degrees of freedom, which 
matches electromagnetism, and the Hamiltonian is posi- 
tive throughout the full physical phase space. However, 
if instead the squared constraint 04 is replaced by the 
equivalent constraint 4 = (B 2 — B 2 — b 2 ) that spans the 
same constraint surface, then a different set of results 
holds. In this case, additional constraints appear from 
the Poisson-bracket relations that are not equivalent to 
the set defined above, and Dirac's counting algorithm 
fails to determine the correct number of degrees of free- 
dom. The resulting theory with 04 replacing 04 is not 
equivalent to the Lagrangian approach. 

Evidently, care must be used in working with a squared 
constraint equation. The constraints 4 and 04 are 
redundant, and the Hamiltonian system is irregular. 
Nonetheless, with these caveats, the KS model with a 
squared Lagrange-multiplier potential provides a useful 
model of spontaneous Lorentz violation. It allows an im- 
plementation of the symmetry breaking that does not 
require enlarging the phase space to include a massive 
mode or nonlinear couplings with A. The only physical 
degrees of freedom in the theory are the NG modes that 
behave as photons. 
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D. Bumblebee Models with (r 2 + r 3 ) / 

In this section, the constraint analysis is applied to 
bumblebee models in flat spacctime that have a La- 
grangian ([7|) with a generalized kinetic term obeying 
(t2 + T3) / 0. Such models do not have a Maxwell form 
for the kinetic term. Throughout this section, arbitrary 
values of n, t 2 , and T3 are used; however, it is assumed 
that discontinuities are avoided when these parameters 
appear in the denominators of equations. The three po- 
tentials in (@| - © are considered, and is assumed to 
be timelike. Since the kinetic term is not of the Maxwell 
form, it is not expected that the NG modes in these types 
of models can be interpreted as photons. For this reason, 
the interaction term B^J^ is omitted in this section. 

The point of view here is that the generalized bum- 
blebee models originate from a vector-tensor theory of 
gravity with spontaneous Lorentz violation induced by 
the potential V. In this context, the vector fields B^ 
have no matter couplings and reduce to sterile fields in 
a flat-spacetime limit. Nonetheless, NG modes and mas- 
sive modes can appear in this limit. Dirac's Hamilto- 
nian analysis is used to examine the constraint structure 
and the number of physical degrees of freedom associated 
with these modes. Comparisons can then be made with 
the results in electromagnetism and the KS bumblebee 
models. 



1. Linear Lagrange- Multiplier Potential 



Beginning with a model with the linear Lagrange- 
multiplier potential in Eq. ((4]), the Lagrangian is given 
in terms of the five fields Bo, Bj, and A. From this the 
conjugate momenta are found to be 



n° = (r 2 + T 3 )(d B Q ) - isidjBj), 

IP = (T 1 -T 2 )(d B j )-T 1 (d j B Q ), 

n< A > = 0. 

The canonical Hamiltonian is then given as 



(47) 
(48) 
(49) 



H = 



- (n - r 2 f 
2(n - r 2 ) 

n 



1 



2(n - t 2 ) 



i W(d j B ) + ~( Tl -T 2 )(d j B k y 



1 



T3 



U°diBi 



2(r 2 +r 3 ) 

T 2 T 3 



(n ) 2 



T 2 +T 3 J - " \2(t 2 +T 3 ) 

+\{Bl-Bj-b 2 ). 



i.djB 3 f 



(50) 



Four constraints are found for this model: 

<h = n< A > 



-(B 2 -B*-b 2 ) 



B„ 



1 



n - t 2 



n 



n - t 2 



-(djBo) 



+B 



T 2 + T 3 



-n u 



T3 



T 2 + T 3 



(51) 
(52) 



(53) 



-X(B 



T 2 + T 3 
T\ - T 2 



(Bj 



T1T3 



2(n - r 2 ) 



tI{t 2 +t 3 ) 

2(n - r 2 )2 

T\T 3 



(d,B ) 2 



\{r 2 



T 2 + T 3 Tl ~ T 2 

T S )Bjd k d k Bj 



(djBj 



- (n + t 3 ){t 2 + r 3 )) 



(ri 



n - t 2 

-r 2 f 



BjdndkBk 



2(n - r 2 ) 

T3 



B djdjB 



2(ri - r 2 ; 
1 



BjidjU") 
1 



~T 2 ) 
T3 



;T 3 



2(n - r 2 ) 
n (r 2 + r 3 



B (djW) 



Tl - T 2 



WO.Bq 



Tl 



Tl 



1 



2(r 2 + t 3 



t 3 2(ri - t 2 ) 

■ (n ) 2 - 



n°{d 3 B 3 



T 3 



2(n - T2 y 



(54) 



The constraint (j)i is primary, while <f> 2l <p 3l and ^>4 are sec- 
ondary. All four are second-class. According to Dirac's 
counting argument there arc n — ni — n 2 /2 = 5 — — 4/2 = 
3 degrees of freedom in this model. 

The constraint <fi 2 shows that only three of the four 
fields Bfj, are independent. In the timelike case, it is 
natural to solve for Bq in terms of Bj. The first and 
third constraints can be used, respectively, to fix n^ A ) to 
zero and to determine II in terms of Bj and IP . The 
remaining constraint ^4 can be used to determine A in 
terms of Bj and IP. Interestingly, this leaves the same 
number of independent degrees of freedom as in the KS 
bumblebee model with a similar potential. One might 
have thought that switching from a Maxwell kinetic term, 
which results in the removal of a primary constraint n° = 
0, would have introduced an additional degree of freedom. 
However, instead, new secondary constraints appear that 
still constrain 11°, though not to zero. As a result, Bq 
and n° remain unphysical degrees of freedom despite the 
change in the kinetic term. 

Since the generalized bumblebee model is not viewed 
as a modified theory of electromagnetism (e.g., no cur- 
rent J M is introduced), there is no analogue or modified 
version of Gauss' law as there is in the KS bumblebee 
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model. Nonetheless, in the constraint ^4, A plays a simi- 
lar role as a nonlinear source term for the other fields as it 
does in the KS bumblebee. Indeed, the constraint equa- 
tion <f>4 w reduces to the same modified form of Gauss' 
law as in (|24|) with J° = in the limit where II — ► 
and the coefficients ti, t 2 , t 3 take Maxwell values. Thus, 
when considering initial values of the independent fields 
Bj and IP in the generalized bumblebee case, the con- 
straint 04 can play a role similar to that of the modified 
Gauss's law in the KS bumblebee model. 

Restrictions on the coefficients r 1; r 2 , t 3 can be found 
by examining the freely propagating modes in the the- 
ory. Investigations along these lines with gravity included 
have been carried out by a number of authors [H, [3(| • 
Since the theory with generalized kinetic terms has three 
degrees of freedom, there can be up to three independent 
propagating modes. These include the NG modes asso- 
ciated with the spontaneous Lorentz breaking. To de- 
termine their behavior, it suffices to work in a linearized 
limit and to look for solutions in the form of harmonic 
waves. Carrying this out in the Hamiltonian formulation 
requires combining the linearized equations of motion to 
form a wave equation for Bj. For physical propagation, 
i.e., to avoid signs in the kinetic term that give rise to 
ghost modes, the condition (n — T2) > must hold [30j ] . 

In this case, two solutions are found that propagate as 
transverse masslcss modes at the speed of light. However, 
a third longitudinal mode can be found as well. In an 
observer frame with wave vector fc M = (ko, 0, 0, £3), it 
obeys a zero-mass dispersion relation of the form 



(Tl - T 2 )k 2 a + (T 2 + T 3 )kl = 0. 



For physical velocities, the ratio 
a = kl/kl 



T 2 + Tj 

Ti - r 2 



(55) 



(56) 



must be positive, which together with the requirement of 
ghost-free propagation gives 



(ti - ti) > 0, (r 2 +r 3 ) < 



(57) 



Note in comparison that the KS bumblebee model has 
(t 2 + r 3 ) = 0, and therefore the third degree of freedom 
does not propagate as a harmonic wave. Instead, it is 
an auxiliary field that mainly affects the static potentials 

0- 

The stability of the theory also depends on whether 
Ti is positive over the full phase space. Examining this 
should include consideration of possible initial values at 
t = that satisfy the constraints. Using integration by 
parts and </> 2 ~ 0, the Hamiltonian (|50|) can be written 
as the sum of two parts, 



tl — Tin -\-TLb- 



(58) 



The first, 



Tin = 



1 



2(71 - T») 
1 



2(r 2 + r 3 ) 



(no + ra^-) 1 



(59) 



includes dependence on the momenta, while the second, 
ti - r 2 



T~Cb = 



(W 2 - 2ri 'I 2 + T3 (^-) 2 



2 w " 2 



(60) 



depends only on the fields B^. 

First consider TCb- From the condition for ghost-free 
propagation in ([57)) , it follows that the first and third 
terms are nonpositive. The second term is nonpositive 
as well if 2ti — r 2 + t 3 > 0, which implies a < 2. Thus 
Tis < if the conditions (|57| hold and a < 2. 

Next consider the momentum-dependent term Tin- 
Assuming the conditions (|57p for ghost-free propagation, 
the first term is nonnegative, while the second is nonpos- 
itive. Note that the two terms are not independent, since 
they are related by constraint 3 . However, one choice of 
initial values that makes both terms vanish (and there- 
fore satisfies 3 « 0) is 



0. 



(61) 



The initial value of A is then chosen to make 04 vanish, 
and B = {b 2 + B 2 ) 1 / 2 is used to make (j> 2 ~ 0. Conse- 
quently, with Tin vanishing, if a < 2, and the condition 
(fST)) holds, then there exist initial conditions with H < 0. 

To investigate the remaining cases, corresponding to 
other possible values of a consistent with (f57|) . use the 
constraint 3 to rewrite Tin as 



1 



Hn = +T 1 d J B ) 2 

2(ti - T 2 ) ( 

[B 3 {W +T 1 d J B )} 2 

— Ol — H 

B 2 



(62) 



In any volume element, choose initial values for Bj of the 
form (B l7 B 2 ,B 3 ) = (0, 0, B{x)). It then follows that 



Tin = 



2(n 



- — - ((n 1 + ndxBof + (u 2 + Tl d 2 B f 
-t 2 ) I 



l- a 



b 2 + B 2 



(n 3 + Tl d 3 B Y 



(63) 



With this form, initial values of the components n 1 and 
n 2 can be chosen that make the first two terms in this 
expression vanish. The third term becomes negative for 
any a > 1, provided an initial value of B 2 is chosen that 
obeys 



B > 



1 



(64) 



With Hn < 0, and n 3 + ti<9 3 5 # 0, the initial value of 
n 3 can then be made arbitrarily large so that the total 
initial Hamiltonian density Ti. = Tin + ris is negative, 
even if TLb > 0. 

Thus, the Hamiltonian density Ti. can take negative 
initial values for any choice of the parameters t± , t 2 , r 3 
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satisfying the conditions (|57| for ghost-free propagation. 
The two examples with a < 2 and a > 1 are sufficient to 
cover all possible cases. 

Evidently a dilemma occurs in the generalized bum- 
blebee model. If the coefficients ti, t 2 , t% are restricted 
to permit ghost-free propagation, then regions of the full 
phase space allowed by the constraints can occur with 
Ti < 0. This parallels the behavior in the KS bumble- 
bee model. With n, r 2 , t 3 equal to Maxwell values, the 
allowed regions of phase space in the KS model include 
solutions with Ti < 0. However, as demonstrated in a 
previous section, if initial values with A = are chosen, 
and current conservation holds, then A = and Ti > 
for all time in the KS bumblebee model. 

Based on this, one could look for similar restrictions of 
the phase space in the case of the generalized bumblebee 
model. For example, the solutions with Ti < described 
above must typically have A ^ at t = to satisfy the 
constraint ^4 « 0. This suggests the idea of trying to 
limit the choice of initial values to A = in an attempt 
to exclude the possibility of solutions with Ti < . 

However, this idea seems unlikely to succeed in the 
case of the generalized bumblebee model, since setting 
A = at t = is not sufficient to restrict the phase space 
to solutions with A = for all time. This is because 
the equation of motion for A has different dependence 
on the other fields in the generalized bumblebee model 
compared to the KS model. In particular, A is not pro- 
portional to just A itself. This is evident even in the 



linearized theory, with B^ expanded as Bp 



Applying the constraint analysis to the linearized theory 
yields a first-order expression for A in terms of £j and W 
equal to 



A 



1 / n + r 3 



26 U1-T2 



(65) 



while the equation of motion for A in the linearized theory 
is 



A 



1 (t"2 + t 3)(ti + r 3 ) 



2b 



(n - r 2 ) 



{d k d k dj£j). 



(66) 



The latter equation shows that (with non-Maxwell values 
t~2 + T3 =/= 0) A is independent of A at linear order. There- 
fore, even if A = at t = 0, nonzero values of A can evolve 
over time. This makes it difficult to decouple regions of 
phase space with Ti > in the generalized bumblebee 
model purely by making a generic choice of initial val- 
ues. It would thus seem likely that the regions of phase 
space with Ti < include solutions obeying A = at 
t = 0. 



2. Quadratic Smooth Potential 

The generalized bumblebee model with a smooth 
quadratic potential |5]| depends on four field components 
Bp and their corresponding conjugate momenta. The 



expressions for n° and IP are the same as in Eqs. (|47j) 
and (pl8"|) , respectively. There are no constraints in this 
model. Thus, according to Dirac's counting algorithm 
there are n~ n\ — ?i2 = 4 — — = 4 independent degrees 
of freedom. This is two more than in electromagnetism, 
and one more than in the KS bumblebee model. 

These four degrees of freedom include three NG modes 
and a massive mode. For arbitrary values of ti, T2, and 
t 3 , all three NG modes can propagate, but with disper- 
sion relations that depend on these coefficients. In con- 
trast, in the KS model, with a Maxwell kinetic term, 
only two of the NG modes propagate as transverse pho- 
tons. A massive mode occurs in either theory when 
V = 2n(B% - B? - b 2 ) 0. In the generalized bum- 
blebee case, there is no analogue of Gauss' law, and it 
is possible for the massive mode to propagate. However, 
in the KS model with a timelike vector, the constraint 
(f3"3")) provides a modified version of Gauss' law, and the 
massive mode is purely an auxiliary field that acts as a 
nonlinear source of charge density in this relation. 

The Hamiltonian for the generalized bumblebee has 
the same form as in (I50|) . but with the potential in the 
last term replaced by the expression in ([5]). With no 
constraints, the full phase space includes solutions with 
an unrestricted range of initial values. Thus, for any 
values of the coefficients t\, t 2 , T3, there will either be 
propagating ghost modes or permissible initial choices for 
the fields and momenta with Ti < 0. 



3. Quadratic Lagrange-Multiplier Potential 

As a final example, the generalized bumblebee model 
with a quadratic Lagrange-multiplier potential ([B]) can 
be considered as well. In this case there are ten fields 
Bp, n M , A, and n^ A ). The conjugate momenta are given 
in |@7]) - (HHJ). The Hamiltonian is the same as in ([BIT]) , 
but with the potential replaced by ©. In this case, two 
constraints are found, 



n (A) 
1 

~2 



(Bl-B)-b 2 f. 



(67) 
(68) 



Constraint 2 imposes the condition @. However, it 
involves a quadratic expression for this condition, and 
therefore the system is irregular, and the same caveats 
must be applied as in the KS model. In particular, sub- 
stitution of an equivalent constraint (f>' 2 — (Bq — B 2 — b 2 ) 
causes Dirac's counting argument to fail. However, with 
4>i and 4>2 identified as first-class constraints, Dirac's al- 
gorithm gives n — n\ — n 2 /2 = 5 — 2 — = 3 degrees of 
freedom. This is again one more than in the KS model. 
In this case, there is no massive mode, and A decouples 
completely. The three independent degrees of freedom 
are the NG modes, which in the generalized bumblebee 
can all propagate. However, even if values of ti, t 2 , and 
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TABLE I: Summary of constraints. Shown for each model are the number of primary (1°), secondary (2°), first-class (FC), and 
second-class (SC) constraints, and the resulting number of independent degrees of freedom (DF). The last column indicates the 
regions of phase space that are ghost-free and have Ti > 0. Current conservation c* M J M = is assumed in the KS models. 



Theory 


Kinetic Term 


Potential V 


Fields 


i° 


2° 


FC 


SC 


DF 


Ghost-Free, H > 


Electromagnetism 




- 




1 


1 


2 





2 


full phase space 


Nambu Model 






A h IF 














3 


subspace (9 3 TP = J°) 


KS Bumblebee 




X(B fl B fl ± b 2 ) 


b^, n M , a, n' A > 


2 


2 





4 


3 


subspace (A = 0) 




(ri = 1, T2 = T3 = 0) 


\n{B^±b 2 ) 2 


b m , n" 


1 


1 





2 


3 


subspace (B^B^ = b 2 ) 






±A(B M B"±6 2 ) 2 


Bp, if 1 , a, n< A > 


2 


2 


2 


2 


2 


full phase space 


General Bumblebee 


non-Maxwell 


\(B„B»±b 2 ) 


Bp, ip, a, n (A) 


1 


3 





4 


3 


no subspace found 




(arbitrary n, ts, T3) 


\n{B^±b 2 ) 2 


b„, n" 














4 


no subspace found 






\\{B^±b 2 ) 2 


Bp, if 1 , a, n (A) 


1 


1 


2 





3 


no subspace found 



T3 can be found that prevent these modes from propagat- 
ing as ghost modes, there are no other constraints in the 
theory that prevent initial-value choices that can yield 
solutions with H < 0. 



IV. SUMMARY & CONCLUSIONS 

Table I summarizes the results of the constraint anal- 
ysis applied to electrodynamics, Nambu's model, the KS 
bumblebee, and the generalized bumblebee. For each of 
the bumblebee models, three types of potentials V are 
considered. The results show that no two models have 
identical constraint structures. In most cases, there are 
one or more additional degrees of freedom in comparison 
to electromagnetism. These extra degrees of freedom are 
important both as possible additional propagating modes 
and in terms of how they alter the initial-value problem. 

In considering the stability of the bumblebee mod- 
els, it is not sufficient to look only at the propagating 
modes. The range of possible initial values must be ex- 
amined as well. In general, when the extra degrees of 
freedom appearing in these models arc allowed access to 
the full phase space, the Hamiltonians are not strictly 
positive definite. However, in the KS models, it is pos- 
sible to choose initial values for the fields and momenta 
that restrict the phase space to ghost-free regions with 
TL > 0. In contrast, in models with generalized kinetic 
terms obeying (t2 + T3) 7^ 0, no such restrictions are 
found. These theories either have propagating ghosts or 
have extra degrees of freedom that evolve in such a way 
that makes it difficult to separate off restricted regions 
of phase space with H > 0. In the end, it appears that 
only the KS models have a simple choice of initial values 
that can yield a physically viable theory in a restricted 
region of phase space. 

The examples considered in this analysis all focused on 
the case of a timelike vector , which is the most widely 
studied case in the literature, since it involves an observer 
frame that maintains rotational invariance. A natural ex- 
tension of this work would be to consider models with a 



spacelike vector . In this case, it is straightforward to 
show that the linearized KS model is equivalent to elec- 
trodynamics in an axial gauge [6j. However, additional 
care is required in conducting a constraint analysis of 
the full nonlinear KS or generalized models, since Bq can 
vanish in the case of a spacelike vector, making addi- 
tional singularities a possibility. Alternatively, an anal- 
ysis in terms of the BRST formalism could be pursued, 
which would be suitable as well for addressing questions 
of quantization. Lastly, an extension of the constraint 
analysis to a curved spacetime in the presence of gravity 
would be relevant, since ultimately bumblebee models are 
of interest not only as effective field theories incorporat- 
ing spontaneous Lorentz violation, but also as modified 
theories of gravity. For example, they are currently one 
of the more widely used models for exploring implica- 
tions of Lorentz violation in gravity and cosmology and 
in seeking alternative explanations of dark matter and 
dark energy. However, performing a constraint analy- 
sis with gravity presents even greater challenges and is 
beyond the scope of this work. 

In summary, the constraint analysis presented here in 
a flat-spacetime limit is useful in seeking insights into 
the nature of theories with spontaneous Lorentz violation 
and what their appropriate interpretations might be. In 
particular, the KS bumblebee models offer the possibility 
that Einstein-Maxwell theory might emerge as a result 
of spontaneous Lorentz breaking instead of through local 
U(l) gauge invariance. Indeed, in the flat-spacetime limit 
of this model, with a timelike vacuum value, electromag- 
netism in a fixed nonlinear gauge is found to emerge in 
a well-defined region of phase space. 
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I 

Erratum Appended to Published Version [Physical Review D 77, 125007 (2008)] 

The Hamiltonian density term TLb given in Eq. (60), in Section III D 1, is incorrect. A correct expression is 

n B = \(n- r 2 ) [(djB k ) 2 - (d,B ) 2 ] - \{r x + r^B. f. (60) 

This term is used to examine the positivity of the Hamiltonian density TL for the case of a model with a general 
kinetic term and a Lagrange- multiplier potential. The change in the term TLb alters some of the conclusions that 
follow Eq. (60), which are stated in terms of a parameter a defined in Eq. (56). A revised argument still assumes 
a > and that Eq. (57) holds for ghost-free propagation. The argument and conclusions for a > 1 are unchanged, 
with the result that TL can be negative. However, for < a < 1, a new examination has to be carried out. Using 
Schwartz inequalities, it is found that TLu > and (djBk) 2 — (djB^) 2 > 0. From the latter condition it follows that 
TLb > — \{ T i + T 3){djBj) 2 . For (ti + r 3 ) > 0, corresponding to a < 1, TLb can be made arbitrarily negative, and 
solutions with TL < can therefore exist. However, for (ti + T3) = 0, corresponding to a = 1, it follows that TLb > 0, 
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and that the Hamiltonian density TL is nonnegative. The model with a = 1 and TL > has a Lagrangian density 
with a kinetic term proportional to C = —^(d^B u )(d ,1 B v ). It has recently been examined in arXiv:0812.1049 by S.M. 
Carroll, T.R. Dulaney, M.I. Gresham, and H. Tarn. It provides a counterexample to the conclusions summarized in 
Table I for the case with general values of t\ , ti , T3 and a linear Lagrange- multiplier potential. The model with the 
same kinetic term but with a quadratic Lagrange-multiplier (as discussed in Section III D 3) also has TL > when 
(ti + T3) = 0. However, all of the other results in the paper, including the discussion of the KS model and the smooth 
quadratic potential, are unchanged by this correction. 



